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摘 要 : 用 有 限 元 通用 程序 进行 稳定 分 析 给 出 的 是 一 个 综合 性 的 计算 结果 ， 探 讨 近似 计算 方法 是 很 有 必要 
的 。 本 文 用 能 量 法 推导 了 在 竖 向 均 布 荷载 作用 下 ， 分 别 考虑 保 向 力 、 非 保 向 力 、 横 向 静 风 荷载 时 
抛物 线 拱 肋 的 稳定 系数 的 理论 计算 公式 ， 并 应 用 于 某 工程 实例 ， 计 算 结果 与 ANSYS 计算 结果 的 
相对 误差 分 别 为 3.4%、2.1%、3.8%。 此 公式 既 简 便 又 符合 实际 情况 ， 可 作为 设计 人 员 对 软件 计 


算 结 果 的 验证 。 
关键 词 : 拱桥 ， 竖 向 均 布 荷载 ， 能 量 法 ， 稳 定 系数 
分 类 号 : AMS(2000) 34D35 中 图 分 类 号 : U443.3 文献 标识 码 : A 
1 引言 


在 以 前 的 研究 中 ， 大 多 只 是 对 圆 弧 拱 受 径 向 荷载 的 情形 进行 分 析 ， 目 的 就 是 为 了 获得 解析 
解 ， 而 对 抛物 线 和 芒 链 线 则 要 进行 试 算 。 贺 弧 拱 在 小 矢 跨 比 情形 下 与 抛物 线 拱 有 较 好 的 近似 
性 ， 所 以 有 学 者 比较 了 圆 弧 拱 、 抛 物 线 及 芯 链 线 拱 的 稳定 有 效 长 度 D,。 而 实际 情况 是 大 多 数 
拱 受 的 荷载 为 竖 向 荷载 ， 与 圆 弧 拱 所 受 的 径 向 荷载 有 较 大 差异 。 本 文 利 用 能 量 法 ， 研 究 在 竖 向 
均 布 荷载 作用 下 ， 把 届 曲 变形 作为 一 个 平衡 状态 来 分 析 并 以 此 来 确定 临界 荷载 8 引 ， 推 导 抛 物 
线 拱 平面 外 横向 稳定 系数 的 实用 计算 方法 。 在 研究 中 分 别 考虑 保 向 力 和 非 保 向 力 的 作用 ， 并 推 
导 了 横向 风 荷 载 作 用 下 拱 的 稳定 计算 公 了 


2 ” 拱 助 的 面 外 稳定 计算 


2.1 抛物线 拱 肋 在 保 向 力作 用 下 的 平面 外 稳定 计算 

所 谓 “ 保 向 力 ”， 就 是 假定 作用 在 拱 轴线 上 的 竖 向 荷载 在 拱 轴 侧 倾 过程 中 始终 保持 紧 直 方 
向 不 变 ， 它 是 为 研究 方便 而 对 荷载 模式 所 作 的 简化 。 

在 保 向 力 荷载 作用 下 平面 外 稳定 研究 过 程 中 ， 我 们 作 如 下 假定 : 

1) 拱 肋 截面 特性 沿 跨 径 方向 不 变 ， 即 等 截面 ; 

2) ” 拱 肋 的 轴 向 应 变 为 零 ; 

3) 作用 在 拱 轴线 上 的 荷载 g 始终 竖 直 向 下 并 沿 跨度 均匀 分 布 。 如 图 1 所 示 。 

设 拱 轴线 的 方程 为 


y = mz, (1) 
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f 为 和 撩 高 ， 工 为 半 跨 长 。 若 记 z = yL, n = Ad. = 2mL, 则 有 


f=mL?, dz= Ldg, (2) 


1 
y= 5nbe. (3) 


8 


Bi 拱 结 构 坐 标 图 


考虑 横向 挠 曲 及 扭转 ， 将 上 面 的 抛物 线 拱 轴 方 程 代入 绕 y 轴 、z 轴 的 挠 曲率 得 
1  dY-6dg nd 


Xv = y ds ds 
= "a + n2g?)-2u, 一 J +n?) hu 十 iU + n292)- 86, (4) 
1 dé Y L- 3 
入 2u! Lo pente 20,, (5) 
RP u, 10, 表示 对 只 的 导数 。 


对 于 承受 对 称 荷 载 的 拱 肋 结构 ， 沿 跨度 水 平方 向 的 位 移 为 零 ， 又 由 拱 肋 的 不 可 伸缩 假定 ， 
拱 的 轴 向 应 变 也 为 零 ， 可 得 


iier (uj (6) 


2mLo 
由 拱 肋 横 向 挠 曲 所 产生 的 应 变 能 为 


xt pda ' El ?LA/1-4 n?g?d (7) 
二 了 yXy Sc m yXy n^" ap, 
由 拱 肋 扭转 所 产生 的 应 变 能 
1 
á GJxĉds = - J GJXZLV1 + n?g? dy, (8) 
一 
荷载 所 做 的 功 为 
L 1 
V= J. wadz=g | vyLdo. (9) 


考虑 到 满足 边界 条 件 ， 取 多 项 式 位 移 函 数 


u=) aip, 6—M bi, i=0,1,2,3,4. (10) 
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对 于 拱 脚 采用 固 结 支 承 ， 边 界 条 件 为 p = X18, u=u, =0, 9 = 0 = 0， 得 
U=a0g, U, = a09%, UU, = 09. (11) 
0 = bog, 6, = bog, 0, = bogs, 


其 中 
g—-1-2z^ +y, gp =-—4p Ap), gi — 412, (12) 


由 式 (6) 可 知 ， 当 yp = 士 1 时 ， 可 以 得 到 
v -483 2 十 3p4 — yf) (13) 
Y 3nL j eed 
通过 将 式 (11) 和 (12) 先 代入 式 (4) 和 (5)， 再 代入 式 (7) 和 (8)， 及 将 式 (13) 代入 式 (9)， 可 
得 总 势能 为 
H-U,4U,-V 
2 
x e LEO 9 (Ap HED) + 95 (As B9) -aiic, (14) 


拱 脚 为 固 结 支承 时 ， 积 分 常数 C1 = 228. 
由 6 了 I = 0， 可 以 得 到 关于 ao 和 bo 的 齐 次 方程 组 ， 则 满足 齐 次 方程 组 有 非 零 解 条 件 的 荷载 即 


为 拱 肋 横向 稳定 临界 荷载 qere 车 记 


EI 
qer = ArT?’ (15) 
AP Ar 称 为 横向 稳定 系数 。 由 此 ， 我 们 可 以 得 到 
| 4n (A2 + kB3)? 
M7 m [Ai tB | (16) 


代入 不 同 的 矢 跨 比 和 刚度 比 ， 由 上 式 即 可 求 得 相应 的 横向 稳定 系数 。 

2.2 抛物线 拱 肋 在 非 保 向 力作 用 下 的 面 外 稳定 计算 

当 拱 肋 侧 倾 时 ， 吊 杆 有 一 个 阻止 侧 倾 的 水 平 恢复 力 妃 。 这 个 恢复 力 就 是 吊 杆 拉力 开 的 水 平 
分 力 ， 它 有 使 拱 回 到 原来 的 平面 内 位 置 的 趋势 ， 这 一 效应 称 为 荷载 的 “ 非 保 向 力 效 应 ”。 由 于 
吊 杆 下 端 受 侧 向 刚度 较 大 的 桥 面 约束 产生 较 显著 的 非 保 向 力 效 应 ， 从 而 使 拱 的 横向 稳定 性 大 大 
提高 四 ， 如 图 2 所 示 。 

拱 肋 在 非 保 向 力 荷 载 作 用 下 ， 平 面 外 稳定 研究 过 程 中 在 2.1 的 基本 假定 基础 上 再 假定 : 桥 面 
系 具有 横向 抗 弯 刚度 EpIp， 和 忽略 其 竖 向 抗 弯 刚度 。 根 据 基本 假定 ， 帅 杆 的 拉力 为 : T= qd. 
其 中 4g 为 桥 面 竖 向 均 布 荷 载 ，d 为 吊 杆 间距 。 水 平 恢复 力 为 


u— up 
H =T 
Ty (17) 
其 中 心 up 分 别 为 拱 肋 和 桥 面 在 z 轴 方 向 的 线 位 移 。 


水 平 恢复 力 所 做 的 功 为 
1 fF 1 q 
Vo og | ecu - [i oy dy E 


第 3 期 谭 红 霞 等 ， 拱 结构 平面 外 稳定 系数 的 实用 计算 方法 499 


图 2:” 非 保 向 力 效 应 


桥 面 系 横向 弯曲 应 变 能 i 
= 3 EpIp(up)? dz, (19) 


AP up = cog 为 桥 面 系 的 横向 位 移 函 数 ，w%p 是 对 z 求 二 阶 导数 。 
由 式 (18) 和 式 (19)， 可 得 总 势能 为 


1 bi 
w = SEI, LS S (A1  KB,) + O(A + kB) + É P (Aa + kB) + 2 kp] 


z [= (ao — co)? Di], (20) 


人 
分 常数 Di = 对 E- 
357 
根据 5IT = 0, 可 以 得 到 临界 荷载 


, Ely 
= eT (2L)?' (21) 


式 中 和; 为 非 保 向 力作 用 下 拱 的 横向 稳定 系数 
A = A B- v B? — 4AC), 


A, B, CHA, rd Aa, Bi m Bs, Di 的 函数 。 

2.3 ”横向 静 风 荷载 作用 下 拱 的 面 外 稳定 计算 

在 考虑 紧 向 荷载 的 基础 上 ， 再 考虑 横向 风 荷 载 的 作用 ， 研 究 拱 肋 的 平面 外 稳定 问题 ， 研 究 
过 程 中 考虑 非 保 向 力 效 应 。 

假定 作用 在 桥 面 系 上 的 横向 静 风 荷载 沿 桥 跨 径 方向 均匀 分 布 ， 其 荷载 集 度 为 gfi; 作用 在 
拱 肋 上 的 横向 风 荷 载 沿 高 度 方向 均匀 分 布 ， 荷 载 集 度 值 qf2 按 拱 肋 0.65 倍 高 度 处 的 风速 值 确 
定 l89]， 风 荷载 垂直 作用 于 拱 肋 与 桥 轴线 所 确定 的 平面 。 

桥 面 系 横向 静 风 荷载 所 做 的 功 


L 1 
Vg = J upgfıdz a upgriL dp = cogri Fi, (22) 
—L —1 
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式 中 
1 
F = Lf (1— 29? + e*)dg. 
一 1 


拱 肋 横 向 风 荷 载 所 做 的 功 
1 
Vr2 = “ant | 0-295) + e) VL+ nàg? dp = aoqfaF5, (23) 
-1 


式 中 
1 
F = Lj (1 — 2? + yt) /1 + n?g? dp. 
ži 


拱 的 总 势能 为 
Ip = I — Vri — Vro 
1 a2 aobo b2 2c? 
= gEl, | T8 (A1 十 Bi) + T3 (42 + kB2) + 7 (A3 + KBs) + TA kpE: | 
-i [aC S (ao € co)? Dı] — cogri i ud 09g p2 F5. (24) 


Hôr = 0， 可 以 得 到 关于 ao, bo, co 的 非 齐 次 方程 组 。 再 由 克 莱 姆 法 则 可 以 求 得 


D, D, D. 
p » £e 25 
Do Ux (29) 


式 中 Do, Da, Ds, D. AA Tr Aa, Bi x Bs, Ci, Di 的 函数 。 
得 到 拱 肋 横向 位 移 与 荷载 的 关系 曲线 


D 
u = agg = a(l — 29? + yt) = p,ü-?e +o"). (26) 


由 于 拱 肋 横向 位 移 的 最 大 值 发 生 在 跨 中 yp = 0 处 ， 故 可 只 考虑 拱 肋 跨 中 截面 横向 位 移 与 荷载 的 
关系 ， 则 跨 中 横向 位 移 uo 可 以 表示 为 
Da 
up = —. (27) 


3 Af 


某 钢 箱 拱桥 矢 跨 比 f/Lo = 1/4.92. WHARE, HRANE ET, = 3.853227 x 
1010N.m2， 抗 扭 刚度 GJ = 2.95105 x 1019 N.m?; 桥 面 为 混凝土 箱 梁 ， 横 向 抗 弯 刚度 EpIp = 
1.515492x101? N.m?。 桥 面 处 的 设计 基准 风速 为 27.483m/s， 风 荷载 集 度 qrı = 1072.445 N/m; 
拱 肋 0.65 倍 高 度 处 的 设计 基准 风速 为 32.314m/s， 风 荷载 集 度 qfa = 1162.498 N/m。 则 拱 
肋 抗 扭 刚度 与 横向 抗 弯 刚 度 比 k = 0.543727， 桥 面 系 横向 抗 弯 刚度 与 拱 肋 横向 抗 弯 刚 度 
tE kp = 39.330462. 

应 用 以 上 所 推导 的 公式 ， 可 以 分 别 计算 出 抛物 线 拱 肋 在 保 向 力作 用 下 、 非 保 向 力作 用 下 和 
横向 风 蓓 载 作 用 下 的 横向 稳定 系数 ， 如 表 1 所 示 。 
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Ais 拱 肋 在 各 种 荷载 作用 下 的 稳定 性 比较 


保 向 力作 用 下 的 Xer 非 保 向 力作 用 下 的 和 er ”横向 风 荷 载 作 用 下 的 拱 顶 横向 位 移 (m) 


本 文 方 法 39.24 126.65 0.132 
ANSYS 37.91 123.99 0.127 
相对 误差 3.496 2.196 3.896 


横向 风 荷 载 作 用 下 拱 顶 横向 位 移 与 竖 向 均 布 荷载 g 的 关系 曲线 如 图 3 左 图 所 示 。 根 据 上 面 的 
计算 结果 ， 当 不 考虑 横向 风力 的 作用 时 ， 其 临界 荷载 为 qu. = 3.29 x 10° (N/m); 当 考 虑 横向 风 
力 的 作用 时 ， 其 临界 荷载 为 ger = 3.13 x 109 (N/m)， 相 对 不 考虑 横向 风力 的 作用 时 略 有 降低 。 

考虑 不 同 的 横向 风力 gr (500 ~ 1500 N/m)， 根 据 公 式 (27) 计算 得 到 的 拱 项 横向 位 移 与 竖 向 
均 布 荷载 4 的 关系 曲线 如 图 3 右 图 所 示 。 由 图 可 以 看 出 ， 随 着 横向 风力 的 gr 的 增加 ， 拱 肋 横 向 
失 稳 的 临界 荷载 减 小 。 


014 qii=1072. 445N/m 
q,,71162. 498N/m 


dom Hb od 位移 四 


拱 顶 横向 位移 四 


0 500000 1000000 1500000 2000000 2500000 3000000 3500000 ] © 500000 1000000 1500000 2000000 2500000 3000000 3500000 
Nm) SN) 


图 3: EH: 实 桥 风 荷载 作用 下 的 荷载 -位 移 曲线 ” 右 图 : 不 同 风 荷 载 作 用 下 的 荷载 -位 移 曲 线 


4 结论 


根据 能 量 原理 推导 出 了 抛物 线 拱 在 竖 向 均 布 荷载 作用 下 ， 不 考虑 非 保 向 力 和 考虑 非 保 向 力 
效应 时 及 在 横向 风 荷 载 作用 下 的 面 外 稳定 的 理论 计算 公式 ， 并 将 其 应 用 于 一 座 实 桥 的 稳定 性 计 
算 ， 其 结果 与 ANSYS 计算 结果 很 接近 ， 这 种 近似 计算 方法 快速 简便 ， 可 作为 工程 人 员 和 设计 
人 员 对 设计 软件 、 设 计 成 果 的 验证 之 用 ， 以 及 作为 工程 加 固 之 前 的 稳定 性 判断 参考 依据 。 
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Practical Computation Method for Out-of-plain Stability Coefficients 
of Arch Structures 


TAN Hong-xia"?, CHEN Zheng-qing!, FENG Zhou-quan! 


(1- Wind Engineering Research Center, Hunan University, Changsha 410082; 
2- College of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105) 


Abstract: The finite element analysis of stability can merely provide a comprehensive result, it is 
necessary to study approximate computation methods. The theoretical formulas are derived by energy 
methods for out-of-plain stability coefficients of parabola arch under the vertical average load, through 
considering the effects of directional loads, non-directional loads and transverse static wind loads, 
respectively. Our results have been applied to an actual project, results show that the relative error 
is respectively 3.4%, 2.196, 3.896 when compared to ANSYS. Our formulas are simple and practical, it 
can be used as a consult with soft for designers. 

Keywords: arch bridge; vertical average load; energy methods; stability coefficient 
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